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Abstract. We briefly review the recently constructed two orbital microscopic model for iron- 
based superconductors based on 5*4 symmetry With this faithful representation of the 
kinematics of the tri-layer FeAs or FeSe structure, the model provides answers and physical 
pictures to fundamental questions related to the robustness of superconductivity and pairing 
symmetry, unifies different families of iron-based superconductors, casts new insight into the 
connections to the other high Tc superconductors, cuprates, and reveals why an s-wave pairing 
can be stabilized by repulsive interactions. Further progresses include that the model provides a 
clean understanding of band reconstruction observed in magnetically ordered states, which is a 
strong support to the kinematics of the Sa model, and captures the essential low energy physics 
of iron-based superconductors based on numerical results from unbiased quantum Monte Carlo 
simulation. 



1. Introduction 

In the past four years, the main theoretical difficulty in understanding iron-based high 
temperature superconductors [21 El [H [5] stems from the fact that the materials are intrinsic 
multi-orbital systems due to the presence of all five 3d-orbitals of iron atoms. Several initial 
attempts to construct simplified models with fewer orbitals was made in an early stage of the 
research^ [71 [8]. While these models could successfully explain some properties, they all had 
serious drawbacks, which resulted in a belief that any microscopic model composed of less than 
all five d-orbitals is insufficient ^ [TO] . 

Recently, very surprisingly, with a proper consideration of the S4 symmetry, the symmetry 
of the building block-the FeAs or FeSe trilayer, we demonstrate that an effective two-orbital 
model is sufficient to capture the kinematics responsible for low energy physics [Ij. The two 
orbital model includes two nearly degenerate and weakly coupled single-orbital parts that can be 
mapped to each other under the ^4 transformation. Not only has the model greatly simplified 
kinematics which essentially can be described by only three important tight-binding hopping 
parameters, but the model also has manifested a clear physical picture to explore many intriguing 
fundamental properties of the materials and has drawn a deep connection to the other high 
temperature superconductors, cuprates. In this paper, we give a quick review of the model by 
emphasizing its unique aspects and summarizing recent numerical and analytic results [TTl IT2]. 



2. Motivation 

Before proceeding to the S4 model, we discuss several main initial motivations which lead to 
the construction of the model. The motivations stem from a few intriguing questions related to 
iron-based superconductors. 

• What is the relation between iron-based superconductors and cuprates? This question has 
undoubtedly been asked by many researchers since iron-pnictides were discovered. In 
the past several years, it has been shown that the two classes of materials share many 
similar properties, including similar layered lattice structures, phase diagrams, strong 
antiferromagnetism in parental compounds, short superconducting coherent lengths and 
so on. However, it is still highly controversial that whether iron-based superconductors 
allow us to finally apply the induction method for high temperature superconductors. 

• How can the two different families of iron-based superconductors, iron-pnictides and iron 
chalcogenides, be unified? Iron-based superconductors include two families, iron- 
pnictides [21 [3l[l] and iron-chalcogenides[5j. They share many intriguing common properties. 
They both have the highest TcS around 50K[13l |3l [HI [151 [IB]- The superconducting gaps 
are close to isotropic around Fermi surfaces[l71 [TB I [T9 l [20 1 1 ^ 122] and the ratio between the 
gap and Tc, 2A/rc, are much larger than the BCS ratio, 3.52, in both families. However, 
the electronic structures in the two families, in particular, the Fermi surface topologies, 
are quite different in the materials reaching the highest Tc. The hole pockets are absent 
in iron-chalcogenides but present in iron-pnictides[17 ( l20 l |2T 1 [22] . The presence of the hole 
pockets has been a necessity for superconductivity in the majority of studies and models 
which deeply depend on the properties of Fermi surfaces [121 [23l IM [25l [M I7f[ [28l l29[ l30] . 
Therefore, the absence of the hole pockets in iron-chalcogenides causes a strong debate over 
whether both families belong to the same category that shares a common superconducting 
mechanism. 

• Why and how can a s-wave pairing state be stabilized if the superconductivity is originated 
from a repulsive interaction? With a d-wave pairing symmetry, such as in cuprates, the 
positive and negative signs of the superconducting order parameter are equally distributed 
in both real and reciprocal space . Therefore, one can argue that the onsite repulsive 
interaction is naturally avoided in a d-wave pairing state. For a s-wave pairing symmetry, 
even in a sign changed 5^ state[24], the argument can not be held in general. While many 
of us use repulsive interactions to obtain a s-wave pairing or 5^, in iron-pnictides, this 
fundamental problem is largely unaddressed theoretically. 

• Why is the superconductivity so robust and how can a single energy scale dominate in pairing 
gap functions in iron-based superconductors? Observed by angle-resolved photoemission 
microscopy (ARPES), a very intriguing property in the SC states of iron-pnictides is that 
the SC gaps around optimal doping on different Fermi surfaces are nearly proportional to 
a simple form factor coskxcosky in reciprocal space. This form factor has been observed 
in both 122[T71 [311 [IHI |32j and lll[33l[33:J families of iron-pnictides. Just like the d-wave 
form factor coskx — cosky in cuprates, such a form factor indicates that the pairing between 
two next nearest neighbour (NNN) iron sites in real space dominates. In a multi-orbital 
model, many theoretical calculations based on weak coupling approaches have shown that 
the gap functions are very sensitive to detailed band structures and vary significantly when 
the doping changes [9l [26l [271 EHl [221 EDj. The robustness of the form factor has been 
argued to favor strong coupling approaches which emphasize electron-electron correlation 
or the effective next nearest neighbour (NNN) antiferromagnetic (AF) exchange coupling 
Jo [351 [36l [371 ESI EHl Sni El [12] as a primary source of the pairing force. However, 
realistically, it is very difficult to imagine there is a single energy scale in local exchange 
interactions if a multi d-orbital model is considered. 



In order to answer the above questions, it is clear that we have to identify microscopic parameters 
that control the essential physics of iron-based superconductors at low energy. 

3. Gauge mapping between the symmetries of superconductivity pairing and 
kinematics 

The idea to address the robustness of possible s-wave pairing in iron-based superconductors 
stems from a paper by Berg, Scalapino and Kivelson[42j, where they showed that the s-wave 
pairing in a two leg-ladder model for iron-pnictides is equivalent to a d-wave pairing in a similar 
model of cuprates. The s-wave and the d-wave can be mapped by a simple gauge transformation. 
Although their idea was carried out in an one-dimensional model, they raised an indispensable 
general point for any dimensional model, that is, the pairing symmetry and kinematics must be 
considered together. In particular, in a system where the short range pairing dominates, the 
phases or the signs of superconducting order parameters can be exchanged with those of the 
local hopping parameters. 

3.1. Gauge Mapping 

A gauge mapping can be generally described as follows. Given a general Hamiltonian including 
a kinetic energy Hq and an interaction energy Hj, we can perform a gauge mapping that changes 
the sign of local Fermionic operators, namely, 

c, ^ i-lf'c, (1) 

where = 0, or 1. It is clear that such a general gauge mapping can be described by an unitary 
operator U. Under the mapping, the hopping terms in can change the sign so that the kinetic 
energy Hamiltonian is changed. However, the eigenvalues of Hq are not changed since it is an 
unitary transformation. It is also easy to see that the interaction Hamiltonian Hj, if it includes 
standard interaction terms, such as local density interaction, spin exchange interaction and so 
on, in general, is not changed under the gauge transformation. Therefore, under the mapping, 

H = Hq + Hj^ = U+HU = ^0 + Hi. (2) 

If the original Hamiltonian H has a superconducting ground state with an superconducting 
order parameter < CiCj >= A^j, the ground state of the new Hamiltonian must have 
< CiCj >= ( — 1)^^+^^ A^j. Therefore, if the pairing is rather local short range, it is possible 
to change order parameter symmetry by choosing a proper gauge mapping. 

3.2. Gauge mapping between d-wave and extended s-wave pairing symmetries in cuprates 
Let's use cuprates as an example to illustrate the point. For cuprates, the kinetic energy Hq is 
dominated by the nearest neighbour(NN) hopping in a square lattice, namely, 

Ho^-tYl 4aCja + ^- (3) 
<ij> 

where < ij > represents the links between two NN sites, and the order parameter of the d-wave 
pairing symmetry in real space is also dominated by coupling two NN sites which is corresponding 
to the d-wave gap form coskx — cosky in reciprocal space, namely, 

< CiCi^^ >= - < CiCi^^ >= Ao (4) 

where x, y are lattice unit vectors along x and y directions respectively. As shown in fig|lj if we 
set a gauge mapping as 





s-wave hopping 
d-wave pairing 



d-wave hopping 
s-wave pairing 



Figure 1. An illustration of the gauge mapping for cuprates. The left part of the figure shows 
the symmetry of the hopping parameters, d-wave pairing symemtry and the simplified Fermi 
surfaces at half filling. The right part shows their corresponding parts after the gauge mapping. 
The gague mapping changes the sign of Fermionic operators at the red sites. 



where ix labels the coordinate along the x direction of the ith site. After the mapping, we have 

(6) 



and 



<ij>y <ij>x 



(7) 



Therefore, the gauge mapping allows us to exchange the symmetries between the hopping 
parameters and the superconducting pairing order parameters. Before the gauge mapping, 
the hopping has a s-wave symmetry and the pairing has a d-wave symmetry. After the gauge 
mapping, the hopping has a d-wave symmetry and the pairing has a s-wave symmetry. It is 
also important to notice that the Fermi surface topology is changed after the gauge mapping 
as shown in fig|l} In a weak coupling approach, the Fermi surface topology has been the key 
ingredient to determine many properties, such as pairing symmetry. If the pairing is driven by 
repulsive interaction or antiferromagnetic fluctuations, the Fermi surface topology before the 
gauge mapping favors a d-wave symmetry while one can argue that the Fermi surface topology 
after the gauge mapping favors a s-wave symmetry. Therefore, the result obtained from a simple 
Fermi surface topology analysis is actually consistent with the gauge mapping result. We will 
revisit this point for iron-based superconductors in the following subsection. 

3.3. Gauge mapping between extended s-wave and d-wave pairing symmetries in iron-based 
superconductors 

The above analysis has been generalized to iron-based superconductors [1 J. A similar gauge 
mapping which maps a s-wave to a d-wave pairing symmetry can be designed as shown in 
figj2]^l, 2, 1\ 2'). We group four neighbour iron sites to form a super site and change the sign of 
Fermionic operators at the red iron sites in figj2]^l'). This transformation maps the Aig 5- wave 
cos{kx)cos{ky) pairing symmetry in the original Fe lattice to a familiar d-wave cosk'^ — cosky 




Figure 2. The illustration of the gauge mapping for iron-based superconductors: (1,2) and 
(1^, 2') sketch hopping parameters and pairing symmetries before and after the gauge mapping 
respectively; (a-h) band structures and Fermi surfaces for two models [36l [10] before and after 
the gauge mapping: (a,e) the Fermi surfaces, (b,f) the band dispersion along the high symmetry 
lines, (c,g) the Fermi surfaces after the gauge transformation, (d,h) the band dispersions along 
the high symmetry lines after the gauge transformation. The hopping parameters can be found 
in the above two references. 

pairing symmetry defined in the two sublattices of the orginal lattice, where {kx^ ky) and (fc^, ky) 
label momentum in Brillouin zones of the origin lattice and sublattice respectively. This is 
actually a straightforward generalization of the transformation used in the study of the one 
dimensional iron ladder model [12]. 

Following the same procedure, we can ask how the kinetic energy changes. In ref.[l], we have 
shown that all the different tight binding models constructed for iron-based superconductors, 
after the gauge mapping, display very striking simplicity: (1) all Fermi surfaces after the gauge 
mapping are relocated around X' ^ the anti- nodal points in a standard d-wave superconducting 
state in the sublattice Brillouin zone; (2) the bands previously located at the different places 
are magically linked in the new gauge setting. In particular, the two bands that contribute 
to electron pockets are nearly degenerate and the bands that contribute to hole pockets are 
remarkably connected to the former ones. 

The first feature seriously challenges the consistency of previous weak coupling analysis 
depending on the difference of Fermi surface topologies for iron-pnictides and iron-chalcogenides. 
After the gauge mapping, based on the new Fermi surface topology, a standard weak coupling 
analysis would naturally favor a robust d-wave superconducting state with respect to the 
sublattice in the presence of repulsive interaction or NN AF coupling[43l Therefore, as 

shown in figjSj if we reversely map to the original gauge setting, the original Hamiltonian must 
have a robust 5-wave pairing symmetry. Namely, based on the new Fermi surface topology, both 



BaFe2As2 




Figure 3. A sketch of a reverse gauge mapping from d-wave to s-wave for both iron-pnictides 
and iron-chalcogenides. 

iron-pnictides and iron-chalcogenides should have 5-wave pairing symmetry before the gauge 
mapping as long as the Fermi pockets are small. 

The second feature is the key feature that leads us to believe that the low energy physics is 
controlled by a two orbital model which is characterized by two nearly degenerate bands. The 
emergence of the universal kinematics in different models suggest that the underlining electronic 
dynamics should be controlled by a much simplified set of parameters. 

4. The construction of a two-orbital model with the ^4 symmetry 

With above observations and the following important facts in the tri-layer FeAs or FeSe structure, 

• The (i-orbitals that form the bands near the Fermi surfaces are strongly hybridized with 
the p-orbitals of As(Se). Since the dx^z ^tnd dyf^ have the largest overlap with the and 
Pyf orbitals, it is natural for us to use dx'z ^tnd dyf^ to construct the model; 

• The local lattice symmetry at an iron site is D2d in which ^4 symmetry is the key symmetry 
operation; 

• The two As(Se) planes are separated in space along c-axis in the tri-layer structure. Since 
there is little coupling between the p orbitals of the two planes and the hopping integrals 
through the p-orbitals are expected to dominate over the direct exchange couplings between 
the d-orbitals themselves, the two orbital model essentially could be decoupled into two 
nearly degenerate one orbital models; 

• The decoupling must be controlled by the ^4 symmetry since ^4 symmetry links two As(Se) 
planes; 

we draw a natural conclusion that the two-orbital model should be classified according to the ^4 
symmetry and each of the two nearly degenerate one-orbital models should have a unit cell with 




Figure 4. A sketch of the ^4 kinematics with dx^z ^tnd dyf^ orbitals. The green balls represent 
As(Se) in the top layer and the grey balls represent those in the bottom layer. The blue and 
red orbitals represent the dx^^ dy'z orbitals. 



two iron sites. A picture of the effective two-orbital model with the 6^4 symmetry is illustrated 
in figjlj The two orbital model is divided into two weakly and coupled one-orbital models. The 
orbital degree of freedom in the two one-orbital models can be called as two components of an 
5^4 iso-spin. The first one includes the dx'z iii the A sublattice and the dy/^ in the B sublattice if 
we use A and B to label two sublattices, and the second one includes the dx'z in the B sublattice 
and the dy/^ in the A sublattice. The first strongly couples to the top As(Se) layer and the 
second couples to the bottom As(Se) layer. We denote q^t to dia to label Fermionic operators 
of the two models that are connected by the ^4 transformation. Then, a general Hamiltonian 
that describes the kinematics is given by, 

^o,S'4 = i^o + Hq,c, (8) 

where includes the two decoupled one-orbital models and i^o.c which also obey the ^4 
symmetry describes the weak coupling between the two models. As shown in[lj, if we keep 
the leading hopping parameters, 

^0 = Yy^^'^sCOskxCOsky - 11] {cl^Cka + dka^ka) 
ka 

+ 2ti{cOskx + COsky){cl^Cka - dka^ka) 

+2tid{coskx - cosky){cl^Cka + dl^dka) 
+4:t2dsinkxSinky{c'l^Ck+Qa - d'^adk+Qa) 
^2t^{cos2kx + cos2ky){cl^Ckcr + d'^^dka) 

+... 

and 

^o,c = '^2tc{coskx + cosky){cl^dkcr + h.c). (10) 
k 



(9) 




Two weakly coupled electronic 
structures: Iso-spin degree of 
freedom linked by S4 synnnnetry! 

Intrinsic doubling of unit cell: An 
electron pocket and hole pocket 
are paired because of the doubling 
of unit cell. 

Doping is near half -filling. 



Hopping parameters: 
t2>ti>IVI, 

t2, t2 has opposite sign. 

The presence of hole pocket is not 
essential: unifying iron-pnictides 
and iron-chalcogenides 



Figure 5. A sketch of the underhning kinematics and band structure of the ^4 model, and a 
hst of the unique and important properties of the model. 



where we limit the hopping parameters up to the third NN (TNN), including NN hoppings, tix, 
tiy, NNN hoppings, ^2, ^2, and TNN hoppings, ^3. The longer range hoppings can be included if 
needed. For convenience, we have defined tu = (^ix + ^iy)/2, hd = (^i^ — ^iy)/2, t2s = (^2 + ^2)72 
and t2d = {t2 - t'^)/2. 

Ignoring tc, each ^4 iso-spin components provides the following two paired band energy 
dispersions due to unit cell doubling. 



/ 2 . 2 • 2 , Ms{coskx^cosky)±tid{coskx-cosky) 
Ee = ek + 4y q^sm^xsin^y + [ — ^ ^ (11) 

/ o . 2 • 9 , Ms{coskx + cosky) ±tid{coskx - cosky) 
Eh^Ck- 4W t^^sm^xsin^y + [ — ^ ^ (12) 

where Ck = 4:t2sCoskxCosky + 2ts{cos2kx + cos2ky) — ji. can capture the electron pocket at 
M points and E^ can capture the hole pocket at F point, ti^, t2s and t2d are the key parameters. 
By just keeping these three parameters, the model is already good enough to capture the main 
characters of the bands contributing to the Fermi surface. Including interactions, the general 



Hamiltonian is given by 



i,a=l,2 i i 

where a — 1^2 labels the ^4 iso-spin, U describes the effective Hubbard repulsion interaction 
within each iso-spin, U' describes the one between them and J'^ describes the effective Hund's 
coupling. Since the two components couple weakly, we may expect U dominates over and 
J^. Then, in the first order approximation, the model could become a single orbital-Hubbard 
model near half-filling. 

In the original paper [Ij, we have suggested that tis breaks symmetry. This statement 
is not absolutely correct because tis and tid can be exchanged by choosing different orbital 
configurations. With the orbital configuration shown in figjlj tis does not break D2d symmetry. 
It is also important to note that the V point and the (tt, tt) point in Brillouin zone here are 
exchangeable because of the doubling of the unit cell. 

5. Numerical results and supporting experimental evidence for the ^4 model 

In figjSj we summarize the important properties in the ^4 model. These properties include 

• For each iso-spin, there is a pair of bands generated by unit cell doubling: one contributes 
a hole pocket and the other contributes an electron pocket. 

• For each ^4 iso-spin, the filling factor is close to half-filling, an important property to allow 
strong spin fluctuation. 

• The unit cell is intrinsically doubled, which is a very different starting point of kinematics 
from previous models. 

• The hole pockets can be suppressed by adjusting tis oi" ^2s- Fixing t2d and ^2^, reducing tis 
leads to flatter dispersion of the band responsible for the hole pockets. Thus, it suppresses 
the hole pockets. Therefore, the model unifies iron-pnictides and iron-chalcogenides. 

• From orbital configuration, we have t2 > \tis\ > 1 and the signs of t2 and ^2 are opposite. 

These characters have important physical consequences. In the following, we discuss some recent 
progresses associated with them. 

5.1. Quantum Monte Carlo (QMC) Simulation on the S4 model 

The microscopic understanding provides a new opportunity to make use of highly controllable 
and unbiased numerical methods to study iron-based superconductors. In a recent paper [llj, we 
have demonstrated that the model in the presence of electron-electron interaction is characterized 
by strong (0, tt) collinear antiferromagentic (C-AFM) fluctuations as observed in iron-based 
superconductors |44 t HS j H6] . Moreover, the pairing with an extended 5- wave symmetry robustly 
dominates over other pairings at low temperature in reasonable parameter region regardless of 
the change of Fermi surface topologies, namely, the presence or absence of hole pockets at F point. 
For both iron-pnictides and iron-chalcogenides, the pairing correlation, the effective pairing 
interaction and the antiferromagnetic correlation strongly increase as the on-site Coulomb 
interaction increases. 

QMC results are summarized in figj6| The competition between the 5-wave and the (i-wave 
in the S'4model is manifested by ti and t2 where ti favors the d-wave while t2 favors the 5-wave. 
A simple phase diagram from QMC simulation is given as the inset in fig. ^h). For iron-based 
superconductors, ti should be roughly equal to ^2/2 and 1^2! is much smaller than ti and ^2- 
Therefore, the 5-wave should be robustly favored. 




Figure 6. (a) Quantum Monte Carlo calculations about the static spin susceptibility for 
the 5^4 model. The inset shows a significant large ti can change the peak from (0,7r) to (tv^tv), 
which signals the competition between NN hopping and NNN hopping, (b) The effective pairing 
interaction as a function of temperature at < n >=1.1, U = 4.0, t2 = 1.2, ^2 = —0.8 for different 
ti. Inset: the competition between Sxy and (1^2 _y2 depends on ti/t2 and —^2/^2 . 



5.2. Oriented gap opening in the C-AFM state: an impact of intrinsic unit cell doubling 
As we have discussed above, one unique feature in the ^4 symmetry construction is that the 
unit cell includes 2 irons. Most of previous studies are based on models with one iron per unit 
cell. In principle, one can argue that models based on both representations should be equivalent. 
However, since the search of a meanfield result depends on the starting point of kinematics, a 
proper starting point is crucial to obtaining a correct meanfield result. 

In a recent paper, we have demonstrated an important impact of the intrinsic unit cell 
doubling[12j. Unlike in a conventional spin density wave (SDW) state where gaps at any two 
points on the Fermi surface that are connected by the ordered SDW wave vector should be 
opened, the band structure in the C-AFM state appears to be reconstructed by the development 
of the magnetic order. ARPES experiments ^47l |48[ |49l [50l |511 152] observe that small pockets or 
hot spots develop around Fermi surfaces. It is fair to say that the scenario of the gap opening due 
to a magnetic order in a multi-orbital system, such as iron-pnictides, is complex. For example, it 
has been argued that the entire Fermi surfaces can be gaped out[53j. Nevertheless, the surprising 
band reconstruction is never fully understood in previous models. 
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Figure 7. The illustration of oriented gap opening in magentically ordered state of iron-based 
superconductors: (a) The C-AFM magnetic state; (b) The illustration of the connections of four 
Fermionic operators in reciprocal space by the C-AFM order; (c) The typical Fermi surfaces in 
the C-AFM state; (d) The band reconstruction and gap opening in the C-AFM state. 



As summarized in figj?} the complicated band reconstruction near Fermi surfaces in the 
magnetically ordered state of iron-pnictides observed by ARPES can be understood in a 
meanfield level by the S4 model. The (0,7r) or (7r,0) C-AFM order does not open gaps between two 
points on theFermi surfaces linked by the ordered wave vector but forces a band reconstruction 
involving four points in unfolded Brillouin zone and gives rise to small pockets or hot spots. 

The doubling of the unit cell leads to an orientation-dependent gap opening in reciprocal 
space. In one T — M direction, the magnetic order splits two folded bands while in the other 
direction the gap can be opened between the pair of electron and hole pocket bands. Such an 
oriented gap opening generates the hot spots or small pockets as observed by ARPES as shown 
in figj7|^c).. 

5.3. Sign change between the top and bottom As(Se) planes in the s-wave state 
The NNN hopping parameters ^2 and ^2 have opposite signs in the ^4 model. This sign difference 
stems from the orbital symmetries. If we take the gauge mapping mentioned earlier, after 
the mapping, the two hopping parameters in the new model have the same signs. Therefore, 
the gauge mapping exactly switches the symmetries between the pairing and the underlining 
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Table 1. A hst of the close connections between iron-based superconductors (iron SCs) and 
cuprates. 

kinematics as shown in figj2]^l, 2, 1\ 2'). 

In the view of the rf-wave picture after the gauge mapping, the superconducting phase for each 
5^4 iso-spin alternates between neighboring squares, which is corresponding to the sign change 
between the top and bottom As/Se planes in the view of the ^4 symmetry. Such an intrinsic sign 
change between the top and bottom As planes along c-axis explains why a s-wave before gauge 
mapping can be stabilized even if the repulsive interaction is responsible for superconductivity. 

The relative phase between the two ^4 iso-spins can take two different values, or tt, in 
principle. The difference can be classified as the irreducible representation of the ^4 symmetry. 
The S/^ point group has four one-dimensional representations, including A, B and 2E. In the 
A state, the ^4 symmetry is maintained. In the B state, the state changes sign under the ^4 
transformation. Therefore, in the A phase, since the ^4 symmetry is not violated, the relative 
phase between the two iso-spin are equal to tt in space, while in the B phase, the relative phase 
is zero. The A phase in the presence of the hole pockets is an sign-changed s-wave or 5^ state. 
Therefore, the model predicts a new superconducting state, B phase state. However, it is still 
an open question whether iron-base superconductors can carry both states or not. 

6. Discussion and perspectives 

6.1. Comparison between cuprates and iron-based superconductors 

By revealing the underlining electronic structure, we can make a close comparison between 
cuprates and iron-based superconductors. In Tablejl} we list the close relations between two high 
Tc superconductors. From the table, it is clear that by determining these physical properties of 
iron-based superconductors listed in the table can help to determine the high Tc superconducting 
mechanism. 



6.2. More supporting evidence for the S4 symmetry 

The S4 model provides a microscopic understanding to two measured universal behaviors: 
the superconducting gap function that is dominated by the single coskxCosky form and the 
almost universal NNN antiferromagnetic exchange coupling strength J2 measured by neutron 
scattering [iS"^ "Ml |55] . The existence of such two single energy scale quantities is also a strong 
support for the ^4 model. If we treat the ^4 model in the presence of strong electron-electron 
correlation, both J2 and the gap function in the s-wave state is determined by the NNN hopping 
parameters. 

The 5^4 model also suggests that the properties of an iron-based superconductor can be 
understood by simply observing bands responsible for electron pockets. Since the physics is 
dominated by the pair of bands that are responsible for a pair of hole and electron pockets. The 



third hole pockets at the F point observed in iron-pnictides is not expected to play a significant 
role. However, the existence of the third hole pocket may change the doping concentration for 
the 5^4 model in a stoichiometric material. Namely, in the view of the 5^4 model, the material 
could be self electron doped. For example, if a large third hole pocket is present, we would 
expect that the 5^4 model could be heavily electron-doped. This expectation naturally explains 
puzzled physical properties of LiFeAs where a large third hole pocket was observed [56]. No 
static magnetic order is observed in the stoichiometric LiFeAs with Tc ^ ISKjbT, ^58]. In the 
view of the 5^4 model, the LiFeAs is self electron doped so that many physical properties of the 
material should be understood from the heavy electron doped side of phase diagram |I59]. 

6.3. Perspectives 

The 5^4 model provides a fresh starting point to explore novel physics in iron-based 
superconductors. Here, we discuss several perspectives for the future research. 

First, since the 5^4 symmetry stems from the fact that As stands out of Fe plane, obvious 
interesting physics may stem from c-axis coupling. Recently, ARPES has provided the existence 
of superconducting gap nodes along c-axis in BaFe2As2-xPx^^^^' It is very likely that the 
5^4 symmetry may be linked to such a development of nodes since the superconducting phase 
along c-axis in the 5^4 model has very intriguing sign change as we discussed earlier. 

Second, electron nematism observed in iron-pnictides [62j can be fully investigated in the 5^4 
model. Namely, the 5^4 model provides a microscopic base to understand electron nematism. 
The 5^4 symmetry breaking can lead to an electron nematic state. In the previous theoretical 
studies, electron nematism is obtained through phenomenological models [63 , IM f [65 } l66j. 

Third, with only a set of three parameters describing kinematics, it becomes possible to 
obtain the dependence of these parameters on lattice parameters. Such an dependence may 
help us to explain empirical observations of the connection between the maximum TcS and the 
As-Fe-As angles[13j. 

Fourth, it is interesting to investigate physical properties on surfaces where the 5^4 symmetry, 
in principle, is naturally broken. A natural consequence of such a symmetry breaking is that 
symmetry broken states may be more stable on surfaces than in bulk. 

Finally, the doping phase diagram and its parameter dependence, as well as many other 
physical quantities need to be investigated in the 5^4 model. 

6.4- Summary 

The 5^4 model provides great simplicity to study and unify the physics of iron-based 
superconductors. It reveals the deep connections between curpates and iron-based 
superconductors. Since the uniqueness of the model stems from the unique structure of FeAs 
or FeSe tri-layer, revealing such an underlining microscopic kinematics cast a new light on the 
search of new structures that can produce high superconducting transition temperatures. 
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